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ABSTRACT 


I 


A study of the unsteady pressure field within the combustor and test section 
of the NASA-Langley Research Center, 8-Foot High Temperature Structures 
Tunnel (HTST) has been performed. Pressure spectra and cross^cor relations-": 
measured at the test section centerline show fluctuations are dominated by tun- 
nel background noise at frequencies below 5 kHz. Above 5 kHz, fluctuations 
display characteristics typical of pressure fields found beneath turbulent bound- 
ary layers. Further analysis of the lower frequency range for hot flow indicates 
a dependence of background noise on combustor conditions. Measurements with- 
in the combustor show a pressure fluctuation with components at 31 Hz and its 
first eight harmonics. The magnitude of this fluctuation seems to correlate well 
with background levels in the test section. For cold flow, ^ the source of pres 
sure fluctuations in the test section was not identified; however, the combustor, 
structural vibration, and the nozzle boundary layer were eliminated as possible 
sources. 


In addition to the main study, the effectiveness of the model protection sys- 
tem in reducing startup and shutdown loads was assessed. Measurements show 
that loads are reduced to values close to those encountered by test articles 
while in the test stream. 
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INTRODUC TION 


The first flights of most airplanes and aerospace vehicles have been pre- 
ceded by development programs which included many hours of scale model wind 
txinnel testing. To apply the information obtained from such tests to full-scale 
vehicles, however, there must be some assurance that a correct reproduction 
of actual flight conditions has been achieved. Unsteady flow, which exists to 
some extent in all wind tunnels, is a possible source of error in such testing. 

By just its influence on boundary layer transition, unsteady flow can affect mea- 
surements of steady lift and drag as well as heat transfer. 


With such possibilities it is advisable to assess the flow quality within a 
wind txannel in order that sources of measurement error can be identified and, 
if possible, eliminated. The investigation embodied in this report examines 
the unsteady pressure field within the 8-foot High Temperature Structures Tun- 
nel (HTST) at NASA-Langley Research Center for this purpose. 

Over portions of a vehicle's surface where the flow is attached and shock- 
free in flight, pressure fluctuations can result only from the turbulent boundary 
layer covering it. The presence of other fluctuations during tunnel testing repre- 
sents a possible source of error. This differentiation between boundary layer k-. 
and nonboundary layer sources of pressure fluctuations provides a natural means 
of identifying spurious contributions to the pressure field within a tunnel. Devi- 
ations from the well documented behavior of boundary layer pressures will 
confirm the presence of these contributions while the specific character of the 
deviations should provide information on their source. 


With such an analysis in mind, overall levels, power spectra and cross - 
correlation functions of pressure measured in the test section are presented. 
Using these, the frequency range over which spurious fluctuations dominate is 
bracketed, and some identification of sources is carried out. 


As a sidelight to the main study, the effectiveness of the HTST model pro- 
tection system in reducing startup and shutdown loads was measured. The 
results of this investigation are presented in Appendix A. 


DESCRIPTION OF EXPERIMENT 


Using a longitudinal line of pressure transducers flush mounted in a flat 



plate, the fluctuating pressure field inside the test stream of the HTST was 
measured for a variety of tunnel operating conditions and two model angles of 
attack (0°, -15°). Concurrently with these measurements, conditions within the 
tunnel combustor section were monitored with a single pressure transducer. 
Boundary layer profile information for the flat plate was derived from traverses 
of a pitot pressure /total temperature probe. Fluctuating pressure data were 
recorded on magnetic tape and later analyzed to determine their rms values, 
power spectra, and cross -correlation characteristics. A description of the 
HTST facility, data acquisition and analysis procedures, and the instrumenta- 
tion employed in the test program is provided in the following paragraphs. 


Tunnel Facility 


The experimental measurements took place in the 8 -ft HTST at Langley 
Research Center. The 8-ft HTST is a hypersonic, blowdown tunnel with an 
open jet test section 2.44m (8 ft.) in diameter and a nominal test section Mach 
number of seven. General details of the tionnel are shown in Figure 1. Working 
fluid of the tunnel is methane-air combustion products. Reservoir conditions 
in the combustor range in temperature from 1366 K to 2367 K (2000°F to 3800°F) 
and in pressure from 6. 894 x 10 N/m^ to 2. 758 x 10^ N/m^ (1000 psi to 4000 > 
psi). If desired, reservoir temperatures in the range of atmospheric conditions 
(300 K) may be obtained by operating the tunnel without combustion. Test arti- 
cles are attached to a retractable sled which can be pitched to angles from -20° 
to +20° with respect to the tunnel centerline. Table I contains a listing of the 
conditions at which the tests in this series were performed. The tabular values 
of local velocity (Uj^ ), dynamic pressure (q^^), and Mach number (M^ ) were mea- 
sured above the test sled surface outside the boundary layer. Free stream 
values q^ , M^) may be determined using oblique shock relations. Depend- 

ing on air supply, test conditions, and nozzle cooling requirements, run times 
varied from 30 to 70 seconds. Further information on the performance charac- 
teristics of the 8-ft HTST may be found in Reference 1. 


Data Acquisition Program 

The essential pressure and boundary layer measuring devices were four 
piezoelectric pressure transducers flush mounted in a flat plate and a remotely 
operated pitot pressure /total temperature probe. These instruments were 
mounted on the tunnel's test panel holder which was injected into the test sec- 
tion when steady flow conditions had been established and retracted when data 
recording was complete. General details of the tunnel test section and panel 
holder are shown in Figure 2. Specificrdimensional information on transducer 
locations is shown in Figures 3 and 4. 
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The transducers were flush, nnpunted in a 45. 7cm x 33.0cm x 2. 5cm 
(18” X 3” X 1”) aluminum plate. Vibration dampers provided shock isolation 
between the plate and the panel holder. A photograph of the plate, showing 
the positions of the dampers, is presented in Figure 5. The plate itself was 
also instrumented with an accelerometer and two thermocouples to measure 
the acceleration and temperature environments experienced by the pressure 
transducers. In addition, static pressure on the plate was measured through 
a port 7. 62cm (3”) off the center line at the same axial station as the third 
pressure transducer. 


A boundary layer trip, consisting of a row of 2.39mm (0.094”) diameter 
spheres, 9. 55mm (0. 374”) apart, was located 12. 75cm from the leading edge 
of the panel holder. Reference 1 indicates that this should have produced a 
fully developed turbulent boundary layer at the transducer locations in all runs 
except No. 4. Transitional flow was expected for this one exception. Two- 
dimensional flow over the panel holder was achieved by attaching side plates 
extending both above and below the panel holder's surfaces. 


The boundary layer probe, pictured in Figure 6, was mounted 174.6cm 
from the leading edge of the panel holder. The sensing elements consisted of 
a pitot pressure probe made frbm stainless steel tubing 1.65mm (0. 065”) in 
diameter and a platinum -platinum rhodium 13% total temperature thermocouple. 
The probe required approximately three seconds to complete a vertical cycle 
(upstroke, downstroke) and extended 7.62cm (3”) above the model surface at 
the top of its stroke. Boundary layer parameters listed in Table I occurred at 
the axial station of the probe. Extrapolation of these values upstream to the 
first transducer position was accomplished using the methods described in 
Reference 2. 


In addition to the instrumentation in the test section, a single pressure 
transducer was placed in the tunnel combustor, upstream of the fuel spray bar, 
to record the level of pressure fluctuations. The general location of this trans- 
ducer is indicated in Figure 7. 


Data Acquisition Instrumentation 

The outputs of all test section transducers were recorded on magnetic tape 
for later analysis. Accelerometer and pressure transducer outputs were am- 
plified prior to recording. Thermocouple outputs as well as boimdary layer 
measurements were digitized at twenty samples /second and recorded in the 
tunnel control room. A block diagram of the data acquisition equipment is 
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shown in Figure 8. Major system components are described in the following 
paragraphs. 


Pressure Transducer (PGB 112M44). The pressure transducers employed 
in the test panel were quartz crystal, piezoelectric devices having a nominal 
sensitivity of 50 mV/psi. Their frequency response was ^3 dB from 10 Hz to 
100 kHz and the diameter of their active sensing area 3.81mm (0. 150"). In the 
present application, the transducers were mounted in water-cooled adapters 
(PCB 64A02) and recessed 0. 127mm (0.005") below the test panel surface. The 
recess was filled with General Electric RTV 560 silicone rubber to reduce heat 
transfer to the crystal through the diaphragm and to decrease any pyroelectric 
effects due to unsteady temperature variations in the tunnel. The thin layer of 
rubber had no noticeable effect upon the frequency response of the transducers. 


Scaling arguments based on the results of previous boundary layer investi- 
gations (Refs. 3-8) indicate the probable presence of pressure fluctuations at 
frequencies as high as 300 kHz within the boundary layer developed in the 8 -ft 
HTST. There are, unfortunately, no commercially available transducers of 
sufficient sensitivity with response this broad. A previous NASA study (Ref. 3) 
conducted at conditions similar to those encountered in the 8 -ft HTST, employed 
a transducer developed specially for that applicatipn. The response of this 
transducer extended from 1 kHz to 300 kHz. However, subsequent analysis , 
showed that the energy measured above 100 kHz had negligible effects upon 
calculated overall pressure levels and boundary layer characteristics. For 
this reason the response of the PCB 112M44 transducer was regarded as ade- 
quate for this investigation. 


Pressure Transducer (PGB 113A24). Measurements of fluctuating pre- 
sure within the combustor were ol^tained with this quartz crystal transducer. 
Its full-scale range was 6.89 x 10 N/m^ (1000 psi) at a nominal sensitivity of 
5 mV/psi. Crystal resonance is specified by the manufacturer as in excess of 
500 kHz. 


Accelerometer (PCB 302A). The acceleration level of the plate was mea- 
sured with this unit. This is a-quartz, piezoelectric de.vice, with a flat response 
from'.5-5000 kHz',- and is capable of operation at temperatures to 395 K: . 

Battery Power Supply (PCB 480A). Each power supply provided a constant 
18 volt, two milliamp excitation current to a single pressure transducer or 
accelerometer. They were used in preference to line power supplies because 
of their low noise characteristics. 
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Thermocouples. Standard chromel-alumel thermocouples measured test 
panel temperatures. They were mounted at the top.'of holes counterboreH to 
within 3. 175mm (0. 125") of the plate surface. Contact between the thermo- 
couples and the aluminiim was maintained by spring pressure. 

Amplifiers (Neff 122). The outputs of the plate pressure transducers and 
the accelerometer were amplified by these units prior to recording onto mag- 
netic tape. They were wideband, differential dc amplifiers which produce 
gains from 0 to 1000 in eleven fixed steps with an additional variable gain from 
1 to 2. 5 at each fixed step. At all values of gain, their frequency response is 
dB from dc to 100 kHz. 

Tape Recorder (Sangamo Sabre 3). ^This is a standard 14-channel instru- 
mentation tape recorder. When supplied with FM wideband group 1 electronics 
and operated at 305 cm/sec (120 ips) it has a flat frequency response from dc to 
80 kHz on all channels. 

Pressure Transducer Calibrator (Photocon PC 125). Prior to each tunnel 
run, this unit was used to check the sensitivity of each pressure transducer in 
the test section. It contains an oscillator, power amplifier, and portable 
driver which will produce a 1000-Hz signal in a continuously adjustable range 
from 100 dB to 160 dB (re 0.00002 N/m ), The frequency of the calibration 
signal may be varied through an external oscillator connection, though at the 
cost of decreased sovind pressure level for lower frequencies. Using this 
connection, sensitivity checks were performed at 250 Hz and at either 130, 140, 
or 150 dB. 


Test Procedure 

Prior to each tunnel run’^ the sensitivity of the pressure transducers was 
checked by placing the portable driver from the acoustic calibrator over each 
unit and recording the resultant voltage due to a known sound pressure level on 
tape. 

Data recording began with pilot burner ignition in the combustor and con- 
tinued until final shutdown of the tunnel after the model was retracted. Thus, 
both the transient pressure loading associated with tunnel starting processes 
as well as the pressure fluctuations beneath the boundary layer during steady 
flow were measured. 


Data Analysis Procedures 

Tape records of pressure and acceleration time histories were analyzed 
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after each r\m in terms of overall levels, power spectra, and cross/correla- 
tions between various pairs of pressure transducers. No correction of the 
pressure records for transducer acceleration sensitivity was required. Spuri- 
ous signal levels generated by acceleration were more than 15 dB below the 
recorded pressure levels at all frequencies. Numerical integration of the digi- 
tized temperature and velocity yielded the boundary layer parameters used to 
scale both spectra and correlation functions. Data were first analyzed with- 
out any filtering of low-frequency data and subsequently reanalyzed with high 
pass filtering at 5 kHz to demonstrate the effect this had on both overall levels 
and the correlation functions. 


Analysis Equipment 

Spectrum Analyzer and Averager (Federal Scientific UA-14 and 1014). 
Power spectra of pressure fluctuations were obtained with this device. Its 
400-point resolution can be applied over frequency ranges from 0-10 Hz to 
0-50 kHz for a real-time auialysis. For recorded signals, the effective anal- , 
ysis range may be increased by playing the tape at some fraction of its record' 
speed. The highest frequency analyzed then becomes a product of the total 
analysis bandwidth set on the machine multiplied by the reciprocal of the tape 
speed fraction. In this manner, the pressure data was analyzed to 100 kHz by 
playing the tape records at half speed and using the 0-50 -kHz range of the 
analyzer. 

Correlator (SAICOR 43A). This instrument accepts dual inputs, A and B, 
and will calculate either the autocorrelation or the cross correlations 

Rab^'*’) ^BA^”^^* Correlations at time delays as small as ZOO nanoseconds 

may be determined. Mean square values of random signals can be determined 
by computing their autocorrelation function and calculating its value at t = 0. 
This method was used to determine overall pressure levels and accelerations. 

Filters (Krohn-Hite 3202). These units are solid-state high- and low-pass 
filters with continuously adjustable cutoff frequencies from 20 Hz to 2 MHz. In 
their maximally flat response mode, they are -3 dB at the cutoff frequency and 
provide 24 dB /octave attenuation outside the pass -band. They were used to fil- 
ter the pressure records prior to the second analysis of this data. 


RESULTS AND DISCUSSION 


During testing and subsequent data analysis, primary importance was 
placed on describing the unsteady pressure field within the test section of the 
HTST, However, the demonstrated effect which combustor pressure had upon 
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test section conditions naturally led to an additional study of the combustor 
environment. The analysis in the following sections parallels this emphasis, 
consideration being given primarily to test section conditions with an analysis 
of the combustor included as an aid to understanding these conditions. 


Test Section 

Several mechanisms exist that could cause pressure fluctuations on the 
surface of an object in the test section of the HTST. Among these are acoustic 
radiation from the combustor and nozzle boundary layer, flow oscillation due to 
total pressure changes upstream of the nozzle throat, and turbulent boimdary 
layers covering the object. Because of the Mach number of the nozzle exit flow, 
there is no possibility that turbulence generated noise (jet noise) radiating from 
the jet mixing region can reach the tunnel centerline within the test section. 

From the viewpoint of reducing measurement errors, it would be desirable 
to eliminate all pressure fluctuation mechanisms except the turbulent boundary 
layer. Since a boundary layer will be present on vehicles in flight, it is the only 
"natural" mechanism and should be present during testing. It is this division 
between boundary-layer and nonboiandary- layer mechanisms which provides 
criteria for analyzing the data in these tests. Deviations of the pressure field 
behavior from that of a turbulent boxmdary layer will confirm the presence of 
additional sources of fluctuation. Further analysis of the specific character of 
these deviations should then give some indication of which source was their cause. 

The behaviour of the pressure field beneath turbulent boundary layers is 
well documented (Refs. 3-8). The random convecting pattern of turbulent eddies 
which causes the velocity perturbations within the boundary layer carries with it 
an associated random pressure field. As the identity of an eddy decays under 
the influence of shear, viscosity, and diffussion, so does its associated pressure 
field. The random nature of the field necessitates that itsudescription be of a 
statistical nature. Explanations of the statistical measures frequencly used in 
this description are given in Appendix B, along with results from past boundary 
layer investigations (Refs. 3-8) which demonstrate their physical interpretation. 
Briefly, the root-mean-square (rms) value is a convenient measure of the aver- 
age dynamic amplitude of the pressure fluctuations, while the power spectrum 
specifies their frequency content and the space -time -correlation f\xnction their 
temporal and spatial coherency. A more rigorous discussion of these parameters 
is found in Reference 9. 

t 

Power spectra, rms values, and correlation functions were computed from 
the pressure data measured in the test section. A comparison between this data 
and results of previous boundary layer investigations (cited in Appendix B) dis- 
plays substantial differences, indicating the presence of a source of pressure 
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fluctuations other than the boundary layer. 


RMS Values. In Figure 9, the rms pressures measured on the plate are 
plotted along with the results from previous investigators. Except in compari- 
son to the results of Chen, the pressure levels for this test range from two to 
five times those measured in other experiments within the same Mach number 
range. However, due to the large scatter of data from other investigations, no 
decisive importance can be attached to this discrepancy. 


Power Spectra. In Figure 10, power spectra of pressure fluctuations 
measured in this experiment are compared to the results of five other investiga- 
tions. The spectra, newly plotted, come from the single cold flow run (No. 11) 
performed and a hot run (No. 12) which displays flow characteristics typical of 
all hot runs (except No. 4, when flow over the transducers was transitional). 
Several areas of disagreement stand out in this comparison. Previous investiga- 
tors have established the spectrum of boundary layer pressure fluctuations to be 
broadband, as one would expect for a quantity whose time history \s random. 


V . j. 2 ' . 

Here, the measured spectra are broadband at high frequencies (rf — > .4), but 

III 


contain discrete frequency components and narrow bands of noise whose ampli- 
tude rises considerably above the broadband level at low frequencies. In the 
cold flow case, four narrow bands of noise occur, centered at 56, 420, 1250, and 
3600 )Hz. In each case, the amplitudes within these bands rise to values twenty 
to fifty times larger than the broadband level expected in this region. For the 
hot flow case, two of the f^our bands of noise reappear (though they are shifted on 
' • ■ ' ■ ’ 6 * 

. the plot because of the different values of in the two runs) and at very low 

frequencies the spectrum is dominated by peaks at 31 Hz and its first eight 
harmonics. This divergence from broadband behavior indicates that the pres- , 
sure field is not random in either hot or cold flow. 


While the bands of noise and discrete frequencies occur at low frequencies, 
the high-frequency portions of the spectra also show discrepancies with previous 
work. The shape of the spectrum appears generally correct in this region (with 
the characteristic high-frequency roll-off) but amplitudes are considerably lower 
than those measured by other investigators. However, this feature can be traced 
directly to the previous discussion of the rms values and the particular way in 
which the spectral curves are plotted. Note that the power spectra G(uj) are 
rnormalized by the ratip/P^^^/.q|- • just the square of the values 

plotted in Figure 8, where they were noted to be slightly high compared to pre- 
vious results. Normalization by this higher ratio leads to low spectral values 

and-the suspicion arises again that the .values P are incorr.ec_tly high. 

tot 1 
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Certainly if their values were reduced, the magnitude of the high-frequency 
components in the spectrum would fall closer to those found in previous inves- 
tigations . 

Cross-Correlation Function. Cross -correlation functions measured her 
tween the leading three transducers are displayed in Figures 11 and 12 for 
both the cold and hot flow runs. Their most striking characteristic is their 
periodicity, a feature absent from results of previous investigations (see 
Figures B-3, 4). Instead of peaking at a single time delay corresponding to 
eddy travel time between transducers, these cross correlations peak at regu- 
lar intervals and little decay of the peak value , occur s . In this situation, it is 
difficult to define a convection velocity. The periodic nature of these plots 
further confirm that the pressure field repeats itself in time rather than fluctu- 
ating randomly. In the cold flow case the correlation function peaks at 
intervals of approximately 800 ^sec, which corresponds to the period of 1,250 
Hz, the center frequency of one of the narrow bands of noise in the spectra. 

For hot flow, the correlation function repeats itself at intervals corresponding 
to the period of 31 Hz, the fundamental discrete frequency already noted in the 
spectrum for that and all other hot flow runs. 

Space Correlation. Space correlation functions for the representative hot 
and cold runs are plotted in Figure 13. The solid lines are the envelope of 
space correlation curves displayed in Figure B-5 from other boundary layer 
studies . 


While the cold flow correlation function seems to be in reasonable agree- 
ment with previous work, the function measured in hot flow is quite different. 
Instead of falling monotonically to a slightly negative value, and then going 
asymptotically to zero as transducer separation increases, the curve stays 
positive over the entire measurement range and performs a small oscillation. 
Because no measurements of space correlation were made at greater transducer 
separation, it is impossible to say whether the correlation function goes to 
zero or remains positive for larger distances. However, even without this in- 
formation, it is apparent that the correlation length defined by this curve must 
be considerably larger with respect to 6* than those measured im:previous work. 
There L <=« 36*. whereas here L equals at least 8&*. If there is little or no 
change in the correlation scale in the transverse direction, the correlation 
area for the hot flow case is about two and one -half times larger than expected. 

Evaluation of Results. The magnitude, spectra, and correlation functions 
computed from the pressure data all indicate substantial deviations from the 
characteristics of a turbulent pressure field. The pressure magnitudes are 
too high, the power spectra are dominated by discrete frequencies or narrow- 
band noise, and the correlation functions indicate that the pressure field no 
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longer has a random character in time. ‘ The pressure correlation areas are 
probably too large, and it is difficult to define either the convection or decay of 
the pressure field. 

I 


The key to identifying the cause of these deviations lies in the power spec- 

/ (iu 4\ t 

tra. They divide into two parts: a low-frequency region (rj — ~ < .4) dominated 

t- 


by narrow-band noise and discrete frequency components, and a high-frequency 
region with the broadband character and roll-off associated with a turbulent 
boundary layer spectrum. It seems probable, that the two parts of the spectrum 
are due to different sources of fluctuation. To test the hypothesis that the high 
frequency components w.ere indeed due to a turbulent boundary layer and that 
their characteristics were being masked beneath the low-frequency fluctuations, 
the data was reanalyzed after being high-pass filtered at 5 kHz to remove the 
lower portions of the spectrum. A frequency of 5 kHz was chosen so that the 
highest frequency narrow band of noise (at 3600 Hz) would be attenuated below 
the broadband level of the remaining portion of the spectrum. 


The results of this second analysis are presented in Figures 14, 15, 16, 
and 17. Comparison with data from previous boundary-layer investigations, 
(Figures B-1 to B-5), now shows substantial agreement. Filtered rms pres- 
sure magnitudes occupy a central position within the scatter of data at those 
Mach numbers, power spectra have not only the approximate shape but also 
the magnitude found in previous work, cross-correlation functions display the 
characteristics of a random pattern with well defined convection speed and decay 
and space correlation functions now show a pattern of pressure which is co- 
herent over nondimensional distances 5/6* compatible with bounda'ry layer 
phenomena. 


The general picture, indicated by the two analyses of the pressure data, 
confirms that pressure fluctuations at frequencies higher than 5 kHz are domin- 
ated by the turbulent boundary on the plate while those fluctuations at lower 
frequencies are dominated by other sources of pressure fluctuation in the tun- 
nel. Insight into the behavior of these other sources can be gained from further 
analysis of the lower.jfrequency region. 


In all probability, several different mechanisms contribute to the lower 
frequency fluctuations. The influence of the boundary layer extends into this 
region, but is obviously masked due to the large magnitudes of the other sources 
There is no particular frequency in the spectra at which a sharp dividing line 
can be drawn between the boundary layer fluctuations and the other sources. 
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Instead, a gradual transition occurs from a frequency region where the bound- 
ary layer dominates to a region where the unknown sources dominate. 
Frequencies near 5 kHz represent the lowest border of the region where the 
boundary layer is the dominant source of fluctuation. 


With this in mind, the total fluctuation of pressure P in the test section 

tot 

will be divided into two parts: Pj^^. composed of the boundary layer dominated 
fluctuations above 5 kHz, and P^, composed of fluctuations from all sources 
below 5 kHz. In Table II, rms values of P^^ are given for most 

of the runs. They represent the average of values measured at the first two 
transducers (data from the third and fourth were sometimes contaminated by 
a crystal resonance problem). 


As Figure 18 demonstrates, P|^j seems to scale on qp a behavior found in 
previous boundary layer studies. The low frequency fluctuationcP^ shows a 
similar behavior in Figure 19 but the scatter about the trend line is quite large. 
A plot of P^ against shows no apparent relation between the two parameters 
and is thus not displayed. Generally, the behavior of P^ shows little depen- 
dence on variation of test section conditions, indicating that its sources are 
probably somewhere else. 


The discrete frequency and narrow band processes which dominate P^^ 
would seem to exclude a variety of aerodynamic sources fr^m consideration 
as the controlling mechanism. Except at very low Reynolc^ numbers, and in 
acoustic feedback situations, aerodynamic noise sources do not display dis- 
crete frequency behavior. Reynolds numbers here are beyond the sensitive 
range. Feedback situations require subsonic flow somewhere to allow acoustic 
excitation to propagate back to the source of disturbances, and at no location 
downstream of the nozzle throat is this condition met except in a thin layer 
close to the nozzle wall. Sound waves propagating upstream within this layer 
will be refracted into the supersonic flow by the severe velocity and tempera- 
ture gradients and then swept downstream. Thus, no feedback loop could 
operate in these conditions. A further argument against an aerodynamic source 
for at least the narrow-band noise lies in the apparent invariance of the band 
center frequencies in substantially different flow conditions. The spectra of 
aerodynamic sources scale on the ratio U/1, where U is a typical velocity, and 
1 a typical length scale in the generating flow. Increasing the flow velocity/^ 
should shift the spectrum to higher frequencies while increasing length scale 
should lower them. From Run 1 1 to Run 12, the velocities increase almost 
threefold, yet no shift in the noise bands occurs. The possible change in 
length scale between the two runs is not enough to cancel this effect. It is 
still possible that the nozzle boundary layer could be radiating the broadband 
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portion of but the dominant discrete frequency components must emanate 
elsewhere. 


Another possible source for the narrow band and discrete frequency fluc- 
tuations which must be eliminated is sound radiation from the structure of the 
nozzle. Complex structures with real edge conditions do not resonate at evenly 
spaced harmonics. Though the structure is undoubtedly excited by the nozzle 
boundary layer, it could not radiate the discrete frequency sound at 31 Hz and 
its harmonics. As far as the narrow band noise goes, more than four natural 
modes of vibration exist for the nozzle below 4000 Hz, and it seems unlikely •; 
that the uniform excitation by the boundary layer in this frequency range should 
cause response only in these four modes. 


Pressure spectra measured in the combustor show strong fluctuations at 
31 Hz and its harmonics, identifying this location as the source of the discrete 
frequency components. Values of rms pressure fluctuation within the com- 
bustor are included in Table II under the heading In Figure 20*, a 

spectrum^f pressure fluctuations within the combustor from 0-300 Hz for 
Tj. 161 L,K ; (2440°F) is compared with the spectrum of fluctuations measured 
during the same run at the first transducer in the panel holder. Except for the 
peak at 155 Hz, each peak in the test section spectrum corresponds to a peak 
in the combustor spectrum. A similar comparison exists in all other hot flow 
runs. Unfortunately, the combustor spectra reveal little about the origins of 
the narrow band noise. In Figure 21, a test section spectrum for T^. = 1611^ K 
covering 0-5 kHz shows the presence of the narrow band noise at 420, 1250, 
and 3600 Hz, while in Figure 22, a combustor spectrum from the same run 
for 0-500 Hz does not show the presence of any narrow band noise at 420 Hz. 
No combustor spectra from 0-5 kHz were obtained for any hot run so the ab- 
sence of the bands at 1250 and 3600 Hz is conjectural. However, in the cold 
flow run, data recorded in both combustor and test section were analyzed from 
0-5 kHz and the comparison of these two spectra in Figure 23 shows narrow 
band noise in the test section but none in the combustor of corresponding fre- 
"8 


*The vertical scale in Figures 20-24 and 28 corresponds to 10 dB SPL per 
inch. However, this holds true only for the first 20 dB below the peak of each 
spectrum, due to the limited linear dynamic range of the spectrum analyzer. 
Below this the vertical scale is nonlinear for each spectrum. Figures 9 and 
14 have been plotted taking this nonlinearity into account, but no compensation 
has been applied in 20-24 and 28 as only qualitative comparisons are being 
made. 
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quency. The apparent matching of peaks in the two spectra at the lowest fre- 
quencies is a result of the inadequate resolution of the analyzer in the 0-5-kHz 
range. Figure 24 shows this feature in the combustor spectrum to consist of 
two peaks, with the test panel peak at 56 Hz falling between them. 


This evidence does not completely exclude the combustor as the source of 
these bands, however. Based upon the manufacturer's specifications, the com- 
bustor pressure transducer was unable to detect pressure changes smaller than 
0.015 psi, or 134 dB sound pressure level. Fluctuations as small as 129 dB 
are large enough to cause the measured levels of narrow band noise in the test 
section if only attenuation of the sound due to wavefront spreading during ex- 
pansion in the nozzle is considered. Thus, the narrow band noise may be 
present, but at a level below the noise floor of the transducer. 


Since the correspondence between the discrete frequency fluctuations in the 
combustor and test section seems established by Figure 20, it is reasonable to 
Expect that the magnitude of will depend upon Pcomb> magnitude of pres- 
sure flucutations in the combustor. In Figure 25, P^ is plotted against 
for the seven runs where values of both are available. Even with this limited 
information, Pu shows a definite increase as the combustor fluctuations increase. 
The lowest value of P^^ is that measured during cold flow. Its appearance close 
on the trend line established by the hot flow data is probably coincidental. It 
should be remembered that Figure 23 shows little relation between combustor 
and test section fluctuations for cold flow; . 


Combustor 

In hot flow, discrete frequency combustor pressure fluctuations are the 
major source of undesirable unsteady pressure variations in the test section. 

The key to understanding test section behavior thus lies in further study of com- 
bustor oscillations and their relation to tunnel operating conditions. 


Values of rms combustor pressure oscillation appear in Table II. In Fig- 

' ■ ' ' ' - ■ . - ■ ^ ^ . , , , , , p 

ure 26, the ratio of fluctuating to steady combustor pressure, comb, is plotted 

versus P^,. Except for data taken during two runs (Nos. 18 and 19), when the 

■ , - - V P 

fuel spray bar cracked, this ratio remains relatively constant at comb 

I’r 

-2 

1.7x10 . The effect of yarying the total temperature within the combustor at 
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constant total pressure upon this ratio'is inconclusive. Figure 27 shows an 


increase of comb with temperature for Pj.~ 6. 8 x 


6 , 2 

10 N/m , 


but a decrease 


with temperature for Pj. ~ 17 x 10 N/m , The limited temperature range en- 
compassed by the data makes conclusions on this subject hazardous. Thus, 
the magnitude of combustor fluctuation appears to depend mainly upon the 
combustor pressure; the fluctuation magnitude increases linearly as the total 
pressure does. 


Combustor fluctuations increase markedly when the fuel spray bar is 
broken. Visual monitoring of the pressure waveform in these conditions also 
shows it to be more impulsive in character than that generated during operation 
with an unbroken spray bar. The power spectra measured during these runs 
display a much greater density of discrete frequency components and ashar:?:’".?' 
monic structure which extends to at least 1000 Hz. A combustor spectrum, 
measured during such conditions (Run 19) is shown in Figure 28. Here, there 
appear to be recognizable peaks every 15-16 Hz across the spectrum. The 
31 -Hz component still has the largest amplitude but energy at approximately 
110, 142, and 248 Hz is quite prominent. 


The appearance of discrete frequency pressure fluctuations within com- 
bustion chambers is a problem often encountered in rocket engine development. 
Reference 10 gives a comprehensive review of the analysis and engineering 
techniques commonly used to solve such problems. The fundamental frequency 
(31 Hz) of the pressure in this case places it in the class of intermediate to 
high-frequency fluctuations as defined in this reference. High-frequency fluc- 
tuations are generally attributed to excitation of acoustic resonances in a 
combustion chamber. Intermediate frequency disturbances do not occur at 
acoustic resonance points but seem to involve coupling of non- re sonant wave 
motion in the chamber with fueMsystem dynamics. (Low-frequency phenomena 
are defined by conditions when the wavelength of interest is much larger than 
typical chamber dimensions. Such is not the case here.) Acoustic resonance 
seems to be the most likely mechanism in this case. Calculation of the funda- 
mental oscillation mode of the combustor,'; taking into account the different 
temperatures upstream and downstream of the spray bar, gives values from 
30. 5 to 33 Hz for the temperature range displayed by Tj. in Table I. The fun- 
damental frequency should vary with combustor temperature, but no consistent 
trend can be found in the combustor spectra from the hot runs because the un- 
certainty in determining this frequency from the spectra is somewhere between 
1 and 1. 5 Hz. Discovery of such a trend would further confirm the acoustic 
resonance hypothesis. 



The spectral structure of the combustor fluctuations during operation with 
a cracked spray bar is not consistent with the acoustic resonance mechanism 
likely during normal conditions. The appearance of energy at 15-16 Hz and 
its odd harmonics is not compatible with the possible standing wave patterns 
within a cylinder closed at both ends. The resonance mechanism does provide 
an explanation of the increased amplitude of the fluctuations, however. The 
amplitude of fluctuations during normal operations are certainly sufficient to 
cause structural vibration in the rings of the spray bar. When a crack appears, 
this flexing of the structure will cause the open area of the cracks to vary peri- 
odically. The methane spilled through the cracks will then also oscillate and 
local combustion processes will follow the supply of fuel. This will feed addi- 
tional energy back into the pressure oscillations and complete the loop. 


CONCLUSIONS 


The unsteady pressure field sensed on the surface of articled placed- in the 
test stream of the 8-ft HTST is separable into two frequency regions within 
which different sources of pressure fluctuations dominate. Above 5 kHz, press 
sure fluctuations are attributable to the turbulent boundary layer on the article; 
below this frequencyj , ^tunnel' background noise is the dominant source. Though 
the bandwidth of the background noise is limited, the total fluctuation energy 
contained in this frequency region is from two to eighteen times larger (de-^'i ’ 
pending upon flow conditions) than the energy contained in the boundary layer 
fluctuations. Thus, the overall character and behavior of the unsteady pres- 
sure field is controlled by the low frequency background noise levels. 


Under hot flow conditions, the majority of interfering noise originates in 
the combustor. Excitation of the fundamental acoustic resonance of the com- 
bustor by the combustion process results in large amplitude (20-40 psi rms) 
discrete frequency pressure fluctuations at approximately 31 Hz and its first 
eight harmonics. These are transmitted to the test section and are dominant 
in the low frequency background noise. 


In both cold and hot flow, the sources of four large amplitude, narrow 
bands of noise have not been identified. Combustion, structural resonance, 
and turbulence generated noise have been eliminated as possibilities. Sound 
radiation from within the air supply system is also an unlikely source, though 
a question exists whether the transducer within the combustor had sufficient 
sensitivity to detect such noise. 





The magnitude of combustion induced pressure fluctuations depends mainly 

upon the total pressure. The ratio of combustor fluctuating pressure to total 

_2 

pressure remained roughly constant at 1. 7 x 10 for total temperatures 
ranging from 1259K (1806°F) to 1889 K ( 2940°F) . 


The model protection system reduces unsteady pressure loads encountered 
during normal startup and shutdown procedures by a factor of approximately 
3. Fluctuating pressure levels are still high, but do not exceed those the 
model encounters while in the tunnel test stream. Details are presented in 
Appendix A. 
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Table II. RMS Values of Fluctuating Pressure Measured in the 
Combustor and Test Section. 
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Figure 1. Overall View, 8 Foot High Temperature Structures Tunnel 
NASA-Langley Research C enter . 












Figure 2. General Configuration of Test Panel Holder 
and Test Section, 8' HTST. 
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1 - Pressure 
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3 - Pressure 

4 - Pressure 
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5 - Accelerometer 

6 - Thermocouple 

7 - Thermocouple 

8 - Static Pressure Port 


Figure 4. Test Panel Detail and 

Instrumentation Locations. 
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Figure 5. Mounting Arrangement of Test Panel. 
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Figure 6. Boundary Layer Probe. 
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Figure 7. Detail of Combustor. 
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Figure 9. Ratio of rms Pressure to Dynamic Pressure Plotted against 
Mach Number. Unfiltered Data from Present Investigation, 
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Figure 11. Unfiltered Pressure Cross Correlations at Stream.’ 
Separations §, Cold Flow (Tj. = 300 K). 
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Unfiltered Pressure Cross Correlations at Streamwise Transduce 
Separations ?»Run #12 (T = 1611 K). 
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Figure 15. Power Spectra of Pressure High Passed at 
5 kHz Compared with Results of Past 
Investigations. 


35 



25 


5 


Pressure Cross Corr 









Figure 17. Filtered Space Correlations Compared with Previous Results 
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Figure 19. 


L/Ow Frequency Fluctuation P Versus Local 
Dynamic Pressure q . 
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Test Section (Vertical Scale Arbitrary), 
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lustor Spectrum, T = 1611 K. 
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Figure 23. Comparison of Combustor and Test Section Spectra for Cold Flow, 
T = 300 K. 
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Figure 24. Comparison of Low Frequency Spectra from Combustor and Test Section 
in Cold Flow, T = 300 K. 
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Figure 26. 


Ratio of Unsteady to Steady Combustor Pressure 
against Steady Combustor Pressure. 
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Figure 27. 


Combustor Pressure Ratio 
Total Temperature. 
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APPENDIX A - MODEL PROTECTION SYSTEM 


During tunnel startup, and prior to establishment of steady flow conditions 
in the test section, the model sled remains in the retracted position shown in 
Figure Al. When test conditions have been established, the model is injected 
into the stream and data are taken. At the end of a run, the model is retracted, 
to keep it out of the collapsing flow in the test section as the combustor and 
ejector are shut off. When it is in the retracted position during the beginning 
and end of a run, the model sled is protected by two panels covered with sound 
absorbing material. This material attenuates unsteady pressure loads on the 
model caused by shock waves and high intensity noise radiation during startup 
and shutdown. The amount of attenuation this material provides has not been 
known in the past. As part of the present test program, a secondary study 
was conducted to measure this attenuation. 


Test conditions in Runs 26 and 29 were repeats of Runs 12 and 14, respec- 
tively. Startup and shutdown procedures were identical within each pair of 
matching runs, except that the protective panels did not cover the retracted 
model infithe. latter run of the pair. Thus, the full fluctuating loads due to • 
startup and shutdown were experienced by the model. These lo^ds were 
recorded by the plate pressure transducers and compared with similar record- 
ings made for the runs when the protective panels were in place. These 
comparisons are shown in Figures A2 and A3, where time histories of the rms 
pressure fluctuations are shown from the beginning of the combustor countdown 
to "model in", and then from "model out" to final stop. In the comparison be- 
tween Runs 12 and 26, it appears that the pressures are attenuated from 8-12 dB 
by the model protection system. In the comparison between Runs 4 and 29, 
attenuations are systematically less and average 6-10 dB. 


Thus, approximately 6 to 12 dB of attenuation is provided by the model 
protection system. The fluctuating loads after attenuation are still high ( 14B+ 
dB), but with fair consistency it seems the protective panels reduce the 
unsteady pressure loads on the model during startup and shutdown to the mag- 
nitude of those it will experience in the test section during a test run. 






Test Chamber 


A-2 
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Figure Al. Position of Acoustic Baffles/Model Protection System 
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APPENDIX B 

STATISTICAL, MEASURES OF A TURBULENT PRESSURE FIELD 


RMS Magnitudes. Because the pressure fluctuations are random, it is 
impossible to predict what their amplitude will be at any specific time. How- 
ever, over long time periods, measures of the average magnitude can be given 
which are called the mean-square value and the root-mean-square (rms) value. 
If p(x, t) is the value of the pressure perturbation at space coordinate x and time 
t, these average quantities are defined as 



mean- square 


root- me an -square 


The rms value is a measure of the average pressure amplitude, while the 
mean-square value is a measure of the average energy of the perturbations. 
Note that both values are strictly defined only in the limit of infinite averaging 
time. In practice, averaging is done over much shorter times with negligible 
resultant errors. 


In Figure Bl, reproduced from Reference 3, a compilation of rms pres- 
sure magnitudes measured beneath turbulent boundary layers by previous 
investigators is displayed. With the exception of one investigation (Chen), 
the ratio of rms pressure to free stream dynamic pressure is generally found 

_ 3 

to be on the order of 10 , with a tendency for this value to decrease with 

increasing Mach number. The two solid lines are semi-empirical predictions 
from Referenced 11 and 12. 


Power Spectra of Pressure Fluctuations. The power spectral density 
function of turbulent pressure describes the general frequency composition of 
the fluctuations in terms of the spectral density of its mean-square value. The 
mean-square value of the pressure in a frequency range between f and (f + AT) 
may be obtained by filtering the pressure with a band-pass filter having sharp 
cutoff characteristics and computing the average of the squared output of the 
filter. This average squared value will approach an exact mean-square value 
as the averaging time approaches infinity. The power spectral density function 
of the pressure measured at location x over the frequency range from f^^ to f 
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is defined by 


^ 2 - 

L.im 


G(f) = Af— fSO ^ p^{x, t, f, Af)df 


It indicates the portion of the total energy in the fluctuations which occurs at 
frequency f. This function is converted to power spectral density at circular 
frequency uj by the relation 


G(uj) 


G(f) 

2TT 


The integral of the power spectral density function is related to the rms mag- 
nitude as 

r 

f. 


rms 


t1/2 


:/ 


G(f)df 


In Figure B2, measured power spectra from several previous investiga- 
tors (References 3 through 8) are plotted. Here the spectra are made nondi- 
mensional '^ith the freestream velocity , the boundary layer displacement 
thickness 6 , and the square of the dynamic pressure q , and then normal- 

2 °° 

ized by the ratio of p /q . The radian frequency uj ( = 2TTf) is made nondimen- 

^00 

sional by the ratio 6/U^ to account for differing boundary layer time scales 
in the various experiments. While the low frequency limits of the measure- 
ments vary, the shapes of the spectra are similar. All have a broadband 
character (no discrete frequencies) extending over two or three decades with 
a sharp roll-off at high frequencies. The low-frequency limits of each data 
set are established by the background noise levels encountered by each in- 
vestigator in the wind tunnel used to perform the experiment. Except for 
Serafini (Reference 5), all investigators have found it advisable to filter out 
frequencies less than 100 Hz. 

Cross-Correlation Function. It is instructive to remember that the turbu- 
lent pressure field is intimately related to the turbulent velocity field within 
the flow. Past hot wire studies in wind tunnels (References 13i 14), boundary 
layers (References 15, 16), and turbulent jets (Reference 17) have indicated 
that the velocity perturbations may be thought of as being caused by a con- 
vecting, decaying, random pattern of eddies. These whorls of fluid, which are 
continuously created by mechanisms of fluid instability, travel with the mean 
flow (though not necessarily at its velocity), gradually losing their identity as 
the' combined influences of shears viscosity, and molecular diffusion extract 
\ and dissipate their energy. Without external inputs of energy into the flow, 
dissipation of the energy in these perturbations results in a decrease in the 
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|fri agnitude of the fluctuations. This is ithe principle of turbulence reducing 
screens used in wind tunnels. A turbulent boundary layer extracts energy from 
the mean flow which borders it. Thus decaying eddies are continuously re- 
placed by new ones created during the process of energy transfer. 


It is reasonable to think of each eddy as having an associated pressure 
perturbation field. As the eddy decays, so must the identity of the pressure 
field. Thus, again there is the picture of a convecting decaying pattern. The 
cross -correlation function of the pressure describes both the spe'ed at which 
the pattern moves, and the rate at which eddies decay and lose their identities. 
The first characteristic is particularly important because a moving pressure 
pattern will cause traveling waves in underlying structures. If the wave speed 
in the structure coincides with the velocity of the pressure field, large struc- 
tural responses occur, with correspondingly increased fatigue and sound 
radiation. 


The cros s -correlation fxinction describes how similar, on the average, the 
pressure measured at one point is to that measured at another point as a func- 
tion of the time delay between observations. Consider two measuring points 
spaced a distance § apart in the streamwise direction. If the eddies in a tur- 
bulent boundary layer were not decaying, so that a frozen pattern of turbulence 
was convecting downstream, the exact pattern observed at the upstream station 
would be observed some time later to pass the downstream station. The dif- 
ference in times at which the same pattern would be observed at the two 
stations corresponds to the time of travel between them. This situation could 
be described by saying that observations at the downstream station are com- 
pletely similar to, or correlated with, those from the upstream station at a 
time delay of t seconds. Here T is the time of travel and is equal to the 
observer spacing § divided by the convection velocity of the pattern U , j 


T 

o 




U 

c 


At time differences greater than or less than T , observations at the down- 

^ o 

stream station are not exactly similar to those registered upstream because 
the upstream observations were of a portion of the pattern which either has 
not arrived or has already passed. 


In the case of decaying, convecting pattern, the concepts of similarity and 
correlation can still be invoked, but it is apparent that observations at two 
stations will never be exactly similar. Some decay in the identity of the pat- 
tern will always occur between the two points. Thus, at the appropriate value 
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of T corresponding to time of travel, there will be maximum similarity Or 
correlation, but not perfect correlation. 


A mathematical function which gives quantitative form to the concepts of 
similarity or correlation between two observations at different times is the 
cross -correlation function. In terms of the pressure p(x, t) measured along 
streamlines, it can be defined 


+ T 


him-’ 

T-^oo 2Tl 


= 


hr S P(-i.t)PU2.t + T)dt 
^ -T 

P (x ) P (x ) 
rms 1 rms Z 


This is a function of the time delay between the signals and its value ranges 
from 1 (perfect correlation) to 0 (no correlation) to - 1 (perfect anti-correla- 
tion). In Figures B3 and B4, cros s -correlation functions measured beneath 
turbulent boundary layers in several previous investigations (References 4, 5, 

7, 8) are displayed. The individual curves shown on each graph represent the 
cross -correlation functions between two points at various separations apart. 
Recalling the picture of a convecting, decaying pattern, a maximum will occur 
in each cross -correlation function at a time delay corresponding to the pat- 
tern's time of travel between them. Because the pattern is decaying, the 
correlation will not be perfect, so that the maximum correlation takes a value 
less than 1. Exactly such features are found in these plots. Each is an orderly 
progression of curves, peaking at successively larger values of time delay 
that correspond to increased spacing, and thus travel time, between observa- 
tion points. Further, a continuing decrease in the value of the peak correla- 
tion occurs with increasing transducer spacing indicating that the eddies decay 
more the further they travel. Note also that each correlation function tends to 
zero rapidly on either side of the peak. This is an indication of the random 
character of the pressure pattern. Those parts of the pattern which convect 
past at time delays less than or greater than the travel time must have little 
or no similarity with parts which cause the maximum or the values of cor- 
relation would remain high. This feature of the correlation functions restates 
the notion that turbulent time histories show no repeatable deterministic 
features either between experiments or within individual measurements. 


The convection speed of the pattern may be determined from these plots 
by dividing the spacing between observers 5 t>y the time delay at which maxi- 
mum correlation occurs. Rate of decay of the pattern eddy may be measured 
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from the decrease of peak correlation with increasing observer separation. 

Space -Correlation Function. Just as the cross -correlation function de- 
scribes the similarity of observations from two points at different times, the 
space -correlation function describes the average similarity of observations 
from two locations at the same time. For two points spaced ^ apart in the 
streamwise direction this function is defined by 


R(§) 


+T 

t'Zc, |p(,= ,,l) P(x2.t)dt 

P (x ) P (x ) 
rms 1 rms Z 


and takes on the same range of values as the cross -correlation function. This 
function provides information both on the effective area over which the pressure 
can be said to act coherently (that is, in phase), and on average phase relation- 
ships between various locations. Two pressures which are perfectly correlated 
must necessarily be in phase, while pressures which are perfectly anti-corre- 
lated must be 180° out of phase. The intermediate value of zero correlation 
indicates a random phase relation between the observation points. 


Knowing the area over which the pressure is approximately in phase is 
important since a structure will respond in different manners if it is being forced 
coherently over a large, rather than small, area of its surface. The way to 
estimate correlation lengths and areas is best understood from examining actur 
ally measured space -correlation functions. In Figure B5, extracted from 
Reference 1, previously measured space-correlation functions of turbulent pres- 
sure are shown for observation points located along lines in the streamwise 
direction. Here the obsjerver separation § has been made nondimens ional by the 
displacement thickness 6 to account for the different scale of the boundary 
layer in each investigation. All of the measured correlations fall rapidly to zero 
within 4 displacement thicknesses. Most then remain slightly negative for long 
distances and are asymptotic to a correlation value of zero. A correlation 
length L can be formally defined by 


L 



R(5)d§ 


but this sometimes results in negligible value of L because of the long negative 
valued "tail" of the correlation function. An alternate, more approximate. 
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measure of L is twice the distance from § = 0 to the point where the correlation 
function crosses zero. By measuring correlation lengths in both the stream- 
wise and transverse direction, their product, called the correlation area, can 
be calculated. This then gives an approximate area over which the pressure 
fluctuates in phase since beyond this region the correlation fxinctions are es- 
sentially zero, indicating a random pahse relationship. 
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Figure Bl. Ratio of rmsiPressure to Dynamic Pressure Plotted 
against Mach Number. Results of Previous Work. 
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